Doping on the kagome lattice: A variational Monte-Carlo study of the t — J model 
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We study doping on the Kagome lattice by exploring the t — J-model with variational Monte- 
Carlo. We use a number of Gutzwiller projected spin-liquid and valence bond-crystal states and 
compare their energies at several system-sizes. We find that introducing mobile holes drives the 
system away from the Spin-Liquid state proposed by Ran et al for the undoped system, towards a 
uniform state with zero-flux. On top of the uniform-state a VBC of the Hastings-type is formed for 
low doping. The results are compared to exact diagonalization on small clusters. This agrees well. 

PACS numbers: 75.10.Kt,75.10.Jm,71.10.Hf 



I. INTRODUCTION 



Spin-liquids were first introduced in 1973 by P. 
W. Anderson as an alternative ground state to the 
anti-ferromagncti. These states were suggested for 
the high- Tc superconductors (such as RVB-states)^. 
Spin-liquids gained tremendous attention due to the dis- 
covery of a number of new materials such as Na4lr308^, 
SrCr8-xGa4-HxOi3ii^, ZnCu3V2(OH)6ClAia and 
KFe3(OH)6(S04)2ii"— , where the realization of such 
states is argued. These compounds share a Kagome like 
lattice in two or three dimensions, with the feature of 
corner sharing triangles (see Fig. [T]), being a prototype 
lattice for frustration with respect for the antiferromag- 
net. The Kagome lattice Heisenberg model has been 
studied extensively and the proposed ground states are 
a valence bond crystali^, a gapped spin-liquid statei^ 
and a particular gapless spin-liquid state (dubbed 
U(l) Dirac state)^'^^. Here the recent study by Yan, 
Huse and White found a gaped Z2 spin-liquid state by 
DMRG, extrapolating remarkable well to the results 
from exact diagonalization,— however a very recent 
rapid communication by Iqbal et al followed such states 
by VMC and found that the U(l) Dirac state, originally 
found by Ran et al still gives the best variational energy 
at half filling. While one may argue the validity of 
VMC due to the bias given by fixing a functional form 
of the wave-function, this can not be argued here, as 
the choice of the wave-function was in fact motivated 
by the DMRG study. One possible reason for this 
puzzling development, might be the difficulty of DMRG 
in two dimensionsi^ii^. We would like to emphasise 
here, that within the VMC the choice of the boundary 
condition changes the variational energy, but not the 
best variational state at half-filling. 
Static impurities on this lattice have been studied by 
Dommange et al^S. This study found the introduction 
of spinless holes leads to dimer-freezing. The interaction 
of the holes was found to be repulsive. A single static or 
mobile hole injected in the Kagome and other frustrated 
lattices has been investigated thereafter with exact 
diagonalization by Lauchli ct al2i finding spin-charge 
separation and dimer-freezing for the positive hopping 



sign of t. Particular the fillings n = 1/3 and n = 2/3 
in the t — J model have been studied with mean-field 
calculations, exact diagonalization and DMRG for a 
number of Kagome like lattices by Indegand et al2^, 
proposing a "bond-order- wave" -state as a generalization 
of a "valence bond solid". In a recent letter by Singh 
the Kagome lattice Heisenberg model with random 
vacancies is studied, and a "valence-bond-glass" state is 
argued^^. The Hubbard model on the Kagome has been 
investigated with dynamic quantum Monte Carlo by 
Bulut et al22 and strong short-range anti-ferromagnetic 
correlation was found by this study. 
In this paper we study mobile holes in the Kagome 
lattice with the t — J-modcl and variational Monte 
Carlo for the first time. The question on doping is of 
very general interest in condensed matter physics. In 
addition very recently doping has been addressed exper- 
imentally on the pyrochlore lattice. This lattice consists 
of alternating Kagome and Triangular planes coupled in 
the direction perpendicular to the planes^. In addition 
the highly degenerate nature of the Kagome lattice, 
makes it interesting to study different perturbation to 
the ground state. This has been pointed out recently by 
one of us in a preprint on the classical model^^ and is a 
viewpoint shared by some of our coUeagues^^. Therefore 
our study has some relevance to the half-filled case too. 
For small thermodynamic doping we test spin-liquid and 
valence-bond crystal trial states previously proposed 
for the Heisenberg model on the Kagome lattice. We 
find that in the presence of a few holes the Dirac 
spin-liquid state is unstable and a zero flux (uniform 
spin liquid) state is formed, accompanied by a "Valence 
Bond Crystal" (VBC). As will be seen below we find 
remarkable well agreement with exact diagonalization 
on small clusters, even we used the most naive boundary 
conditions (periodic boundary conditions). 



II. MODEL AND METHOD 

We study the t — J model within the subspace of single 
occupied lattice sites on the Kagome lattice (see Fig. [T|): 



FIG. 1: The Kagome lattice consisting of corner sharing tri- 
angles. The 3-site physical unit-cell is indicated with dotted 
lines. The underlying Bravais lattice is triangular, leading to 
a hexagonal shaped Brillouin-zone as in the triangular lattice. 
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Here Cja is the electron annihilation operator of an 
electron with spin a on site i, Si is the spin- 1/2 oper- 
ator at site i, and n.i„ ~ cl^Ci^. The sum is over 
the nearest neighbors (n.n.) pairs on the Kagome lattice. 
Vg = Yiii^ ~ "■it"'ji) Gutziller projection opera- 
tor enforcing the single occupied site constraint. We fix 
J = 0.4t in what follows. As wavefunction we consider a 
combination of previously suggested Gutzwiller projected 
MF-states which are classified by the fluxes through the 
triangles and hexagonal respectively, and two types of ex- 
plicit VBCs. Such states are constructed from the mean- 
field hamiltonian: 

Hmf = E X*j4<TCi,<T + h.c. (2) 

This form can be obtained by rewriting the spin- 
operator in fermion operators and introducing the mean- 
field variable Xij = {c\n'^j,iy)- One can define a flux (j) 
through a given plaquette with exp (i(/)) = Opiaquottc Xi] 
with Xij being defined on one edge of the plaquette. 
In addition to the phase (/> on a plaquette, the ampli- 
tude of Xij can be defined on any bond, the latter de- 
fines a VBC. Choosing a mean-field unit-cell capable to 
accommodate the background-flux and VBC, one can 
write out this Hamiltonian for in matrix form. This 
matrix can then be diagonalized numerically to obtain 
the single-particle cigenfunctions. This being the basis 
for our many-particle wave-function in the VMC cal- 
culation, formulated in the usual form as the product 
of two Slater determinants (one for the up- and one 
for the down-spin electrons), where each element is one 
such single-particle eigenfunction. The Gutzwiller pro- 
jection operator defined above is then applied on this 
state leading to the final VMC state to be optimized: 
|*VMc(x»i)) = ^el*MF(Xij))- Projecting is a crucial 



FIG. 2: The flux states U, D and C, being deflned by the 
value of a flux when surounding a plaquette, as explained in 
the text. 

point, as this changes the mean- field energies dramati- 
cally. While this projector can not be treated properly 
in an analytic mean-field treatment, in VMC projection 
is exact. The two VBCs we are choosing arc the one 
suggested by Hastings^^ and the one suggested by Zcng 
and Marston^^. A particular fiux state is parametrized 
with bond-modulation in the pattern of these two VBCs 
(see below), combining background-flux and VBC in one 
state. We use unit cells of the size of 12 and 36 sites, 
corresponding to smaller physical unit-cells of 3 and 9 
sites respectively to accommodate these states. In the 
remainder of the paper we choose the following notation 
for the trial-states: 

1. U: the uniform state with no fiux through the tri- 
angles and hexagons, with a spinon Fermi-surface, 

2. D: the Dirac U(l) state, with zero flux through tri- 
angles and TT-flux through the hexagonals. having 
two Dirac-points instead of a Fermi-surface, 

3. VBC-l(x): the VBC state suggested by Hastings, 
with a 12-site unit-cell, the value x in bracket 
parametrizes the strength of the bond-modulation, 
as the amplitude ratio between strong and weak 
bonds. 

4. VBC-2(x): the VBC patten suggested by Zcng and 
Marston, requiring a 36 site unit cell. The value x 
has the same meaning as above. 

5. C(0): the Chiral state defined by the perturbation 
towards the D-state. Here (/) labels the flux through 
the triangles, while the flux through the hexagons 

is TT — 2(f). 

E.g. U-fVBC-l(l.l) means therefore a state with uni- 
form flux and a bond-order patten resembling the one of 
VBC-1 with a ratio of the strong compared to the weak 
bonds of X — l-l- Similarly we will label other states. 
The states U, D and C are visualized in Fig. [2] the two 
VBCs can be found in Fig. [3l All bonds belonging to the 
patten of the VBC-state in question are therefore given 
the same weight. 

We use a VMC scheme, with the more efficient al- 
gorithm first introduced by Ceperly et al^, this algo- 
rithm computes the ratio of the desired determinants for 
each MC step, with the help of the inverse Matrix. We 
performed a few runs with the naive approach, comput- 
ing directly both determinants to make sure the round- 
ing error resulting from the improved algorithm is not 
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FIG. 3: Left: The patten of the VBC-1 state. The dotted 
hnes indicate the 12-site unit cell. Right: The patten of the 
VBC-2 state. The dotted lines indicate the 36-site unit cell. 
The variational parameter x tunes the ratio between the weak 
and strong bonds. 
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causing any bias. We found that the usual updating 
scheme over n.n. neighbors converges slowly, particu- 
larly at half-filling. The reason is that the weights of 
the wave-function fluctuate strongly, causing the system 
to get trapped in certain configurations for a long time. 
Sampling over n.n. and n.n.n. neighbors improved the 
convergence. We compared the result of both versions 
for special cases and did not find any difference. For 
fractional fillings, it is hard to avoid degenerate states. 
As the result shows very little dependence on the bound- 
ary conditions, we chose to work with periodic bound- 
ary conditions (PBC). This choice is not ideal at e.g. 
half-filling, however comparing different geometries and 
fillings, we found it the most "consistent" choice. One 
sweep is defined as an update over all n.n. and n.n.n. 
neighbors. Typically we performed around 20.000 such 
sweeps to equilibriate the system and another 20.000 for 
measurements. We then average our data over at least 
8 independent runs, and obtained an error-bar. To in- 
vestigate the finite-size behavior we used lattices from 
= 192 to iV = 432 sites. Our larger sizes used, com- 
pare to the sizes currently the larger sizes for VMC on 
the Hcisenberg-modeli^iii. The t — J model is numer- 
ically about one magnitute harder. Because we used 
two physical geometries we will describe the system-size 
by N not by L to avoid confusion as iV = x 3 or 
iV = X 9 depending on the unit-cell in question. Our 
program is flexible in terms of wave- function and lattices, 
this gives us the possibility to check our code against 
published results. To estimate the quality of our result 
we compare our result with the iV = 21 site ED-result 
of Lauchli^. To optimize our variational parameter, we 
used the naive "mesh" approach. We are aware of more 
sophisticated methods such as the SR-method^. The 
later method is undoubtedly very good for the optimiza- 
tion of wave-functions with several variational parame- 
ters, for this flrst exploratory investigation, we chose to 
keep the number of variational parameter small, for this 
reason the mesh- method is sufficient. 



FIG. 4: (Color online) Ground state energy for various states 
and system-sizes in their optimum as a function of doping x 
and the comparison with the ED^ on small clusters of A'^ = 
21. The error-bar is smaller than the symbol-size. 



III. RESULTS 

In Fig. [4] we show the optimized energies of the 
most important trial states function of doping. 

At half filling we recover the result of Ran et ali^ 
the D-state. Already at very low doping (the lowest 
doping-level possible for a spin-balanced state and for 
the sizes investigated) the fiux of this state vanishes, 
as can be seen by comparing the variational energies 
of the U- and D-states. On top of the U-state the 
variational parameters for the VBC-1 and VBC-2 patten 
increase and further lower the energy. The energy 
difference between the pure U-state and the U-|- VBC-1 
and U+VBC-2 states and the variational parameter 
corresponding to the states in questions are plotted in 
Fig. [5l Apparently the VBC-2 patten shows only a very 
small energy-gain compared to the pure U-state and a 
relative big finite-size dependence. The energy difference 
essentially vanishes for N = 324. The VBC-1 patten 
seems to have a stable flnite-size behavior starting from 
N = 300. For higher doping in all cases the VBC-patten 
vanishes. As there is no systematic study on the ground 
state energy as a function of the doping published, we 
included some results on small clusters obtained by 
exact diagonalization by Lauchli ct al (see Fig. [5]). Our 
result is remarkable close to the ED result for N = 21. 
The results from ED do not allow an extrapolation, it is 
therefore not clear how well our energy agrees, yet the 
close match of the two methods which is almost within 2 
percent is very encouraging. At a fllling of 2/3 our result 
shows larger derivation from the ED result, suggesting 
that the real ground state is not described well with our 
wave-function starting from this doping level. 

In Fig. [HI the contributions of the 3 terms in the t — J 
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FIG. 5: (Color online) Above: Difference between the U state 
energy compared to the U+VBC-1 state and U+VBC-2 state 
energies for various sizes. Below: Variational parameter x 
described in the text) for the U+VBC-1 and U+VBC-2 states 
for some sizes as a function of doping x. For clearity we show 
only a few typical error-bars. 



FIG. 6: (Color online) Contributions of the various terms in 
the t — J model (kinetic energy, spin-exchange and "'J'^ -term, 
compare Eq. [T]), for the D+VBC-1 state and the U+VBC-1 
state with the best energy at each doping-level x. 



model for the D+VBC-1 and U+VBC-1 states are plot- 
ted. As can be seen the kinetic energy grows with the 
doping for the U-state and D-state, as we introduce holes. 
The spin-exchange reduces as expected. This reduction 
is relatively small compared for the introduced doping. 
At half filling the D-state has a slightly lower spin-spin 
exchange. At a doping level above x = 0.18 the U-state 
spin-spin exchange energy starts to be lower compared 
to the one of the D-state. The biggest difference between 
the two states comes from the kinetic energy, which de- 
creases much faster for the U-state compared to the D- 
state. This is the reason why only at or extremely close to 
half-filling the D-state is the best variational state. The 
holes remain isolated as can be seen from the nirij term, 
indicating that no phase-separation is taking place. Ap- 
parently the holes move in the background of an almost 
ideal VBC patten for the doping-levels it is formed. The 
low value of the bond-order weight as a variational pa- 
rameter seems to be necessary to keep the holes mobile. 
From the data we can infer the two competing trends 
of hole-mobility and dimerization. Thus in agreement 
with earlier ED studies of static holes and single mobile 
holes, and with the results obtained for 2/3 filling of the 
t— J model. Introducing chirality (C-state) as described 
above increases the energy. The U+VBC-2 state gives 
almost the same energy compared to the pure U-state. 
The D+VBC-2 state gives very similar energies as the 
D+VBC-1 state (which by itself does not differ much 
from the pure D-state). These states are not shown here. 
The variational parameter for the bond-amplitude devel- 
ops a maximum for a doping of around x ~ 0.09 then 
the parameter decreases again. The maximum for the 
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FIG. 7: (Color online) Phase diagram of the t — J model on 
the Kagome lattice. 



1 state is at around x ^ 0.14. Considering the results for 
the energies and the variational parameter, we suggest a 
phasediagram as depicted in Fig. [T] 

To summarize, we have for the first time studied the 
t — J model on the Kagome lattice systematically for 
low doping. Our results are consistent with present ED 
on the = 21 site cluster. We find that the spin-liquid 
state found for half-filling is unstable even for very small 
doping and a zero-flux state is formed. This state shows 
VBC tendencies which vanish for larger doping. There 
are a number of interesting question remaining for future 
studies: most promising are a more free parametrisation 
for the VBCs and the introduction of pairing channels. 
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